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Summary 

A framework is presented for the design and analy- 
sis of quantum mechanical algorithms, the O(jN) step 
quantum search algorithm is an immediate consequence 
of this framework. It leads to several other search-type 
applications - an example is presented where the Walsh- 
Hadamard (W-H) transform of the quantum search algo- 
rithm is replaced by another transform tailored to the 
parameters of the problem. Also, it leads to quantum 
mechanical algorithms for problems not immediately 
connected with search - two such algorithms are pre- 
sented for calculating the mean and median of statistical 
distributions. In order to classically estimate either the 
mean or median of a given distribution to a precision e , 

needs £2(e ) steps. The best known quantum mechan- 
ical algorithm for estimating the median takes 0(e l ) 
steps, and that for estimating the mean takes 0(£ ' ) 

steps. This paper presents 0(E ) step algorithms for 
both problems (all bounds are upto polylogarithmic fac- 
tors). Both algorithms are considerably simpler than 
known algorithms. 

1. Introduction 

1.0 Background Quantum computers were first consid- 
ered in the 70's and early 80's. In 1980, it was shown 
that anything that could be computed by a classical 
computer could, in principle, be computed by a quantum 
computer [Benioff]. But could a quantum computer do 
anything that a classical computer could not? Feynman 
[Feynman] discovered an example of a quantum 
mechanical system which seemed very difficult to simu- 
late on a classical computer, but which a quantum com- 
puter could easily do in polynomial time. 



Through the late 80's and early 90's, the descrip- 
tion of quantum computers was formalized and a num- 
ber of contrived situations and oracle-based examples 
were discovered that a quantum computer could solve 
more efficiently than a classical computer [DJ][Yao] 
[Simon] [BV]. 

The question still remained as to whether a quan- 
tum computer could solve an actual problem of interest 
to computer scientists more efficiently than a classical 
computer. Building a quantum computer is a challeng- 
ing task and before attempting to build one, the design- 
ers would need to be convinced of its ultimate 
usefulness. In 1994, Peter Shor presented an efficient 
quantum mechanical algorithm for factorization - this 
was a problem of great interest to computer scientists 
and it aroused a lot of excitement [Factor]. It was 
expected that this would be quickly followed by other 
quantum mechanical algorithms; however, the next sig- 
nificant result had to wait till 1996 when an efficient 
algorithm for exhaustive search was discovered 
[Searchl]. Since then fast quantum algorithms for other 
important computer science problems have been discov- 
ered [Collision] [Mean] [Median] [Struct], but there is 
still no general technique for obtaining efficient quan- 
tum algorithms. 

1.1 This paper This paper gives a general technique for 
deriving a class of fast quantum mechanical algorithms. 
The idea is to first consider a single quantum mechani- 
cal operation, or a combination of such operations, due 
to which there is a certain probability for the system to 
reach a certain target state t. It is shown that the proba- 
bility of reaching the target state can be made to grow 
quadratically with the number of iterations by iterating 
the quantum mechanical operations in the prescribed 
way. 

The most obvious application of this technique is 
exhaustive search. In the exhaustive search problem, a 
function f(x), x = 1, 2...N , is given which is known 
to be non-zero at a single (unknown) value of x , say Xq 

- the goal is to find x . If there was no other information 

about f(x) and one were using a classical computer, it 

N 

is easy to see that on the average it would take — func- 



tion evaluations to solve this problem successfully. 
However, quantum mechanical systems can explore 
multiple states simultaneously and there is no clear 
lower bound on how fast this could be done. [BBBV] 
showed by using subtle arguments about unitary trans- 
forms that it could not be done in fewer than Q.(Jn) 
steps - subsequently [Search 1] found an algorithm that 

took precisely 0(jN) steps. [Searchl] was derived 
using the W-H transform as the state transition matrix 
and it appeared to be a consequence of the special prop- 
erties of this transform, this paper shows that similar 
results are obtained by substituting any unitary transfor- 
mation in place of the W-H transform. This flexibility 
leads to other search-type algorithms. 

It also leads to applications that, at first sight, are 
not connected to search - two such applications to statis- 
tics are presented. In these a unitary transformation is 
designed due to which the amplitude in a particular state 
comes out to be proportional to the statistic we want to 
estimate. Then by repeating this transformation in the 
prescribed manner, the probability in the desired state is 
increased to a detectable level. Finally a measurement is 
made to determine if the system is indeed in the desired 
state. By repeating the entire operation sequence a few 
times and counting the number of observations, one can 
estimate the probability of occurrence of the desired 
state and hence the original statistic. 

2. Quantum operations 

In a quantum computer, the logic circuitry and time 
steps are classical, only the memory bits that hold the 
variables are in quantum superpositions - these are 
called qubits. Quantum mechanical operations that can 
be carried out in a controlled way are unitary operations 
that act on a small number of qubits in each step. A 
good starting point to think of quantum mechanical 
algorithms is probabilistic algorithms [BV] (e.g. simu- 
lated annealing). In these algorithms, instead of having 
the system in a specified state, it is in a distribution over 
various states with a certain probability of being in each 
state. At each step, there is a certain probability of mak- 
ing a transition from one state to another. The evolution 
of the system is obtained by premultiplying this proba- 
bility vector (that describes the distribution of probabili- 
ties over various states) by a state transition matrix. 
Knowing the initial distribution and the state transition 
matrix, it is possible in principle to calculate the distri- 
bution at any instant in time. 

Just like classical probabilistic algorithms, 
quantum mechanical algorithms work with a probability 
distribution over various states. However, unlike classi- 
cal systems, the probability vector does not completely 
describe the system. In order to completely describe the 



system we need the amplitude in each state which is a 
complex number. The evolution of the system is 
obtained by premultiplying this amplitude vector (that 
describes the distribution of amplitudes over various 
states) by a transition matrix, the entries of which are 
complex in general. The probabilities in any state are 
given by the square of the absolute values of the ampli- 
tude in that state. It can be shown that in order to con- 
serve probabilities, the state transition matrix has to be 
unitary [BV]. The machinery of quantum mechanical 
algorithms is illustrated by discussing two elementary 
unitary operations - these are the W-H transformation 
operation and the selective rotation of the amplitudes of 
certain states. 

A basic operation in quantum computing is the 
operation M performed on a single qubit - this is repre- 
sented by the following unitary matrix: 



M = 



1 



1 



72 

superposition: 



- the state is transformed into the 



J_ J_ 

72 72 



Similarly state 1 is trans- 



formed into the superposition 



1 1 , , 
— , — - . A system 

72 72> 



consisting of n qubits has N = 2 basis states (in this 
paper whenever "a state of the system" is mentioned, it 
will mean one of the N basis states in which each qubit 
is either a or a 1 ). We can perform the transformation 
M on each qubit independently in sequence thus chang- 
ing the state of the system. The state transition matrix 
representing this operation will be of dimension 

2 n x 2 n . Consider a case when the starting state is one 

of the 2" basis states, i.e. a state described by a general 
string of n binary digits composed of some Os and some 
Is. The result of performing the transformation M on 
each qubit will be a superposition of states consisting of 
all possible n bit binary strings with amplitude of each 

n 
2 

state being ±2 . This transformation is referred to as 
the W-H transformation [DJ] and denoted by W . This 
operation (or a closely related operation called the Fou- 
rier Transformation [Factor]) is one of the things that 
makes quantum mechanical algorithms more powerful 
than classical algorithms and forms the basis for most 
significant quantum mechanical algorithms. 

The other transformation that we need is the 
selective rotation of the phase of the amplitude in certain 
states. The transformation matrix describing this for a 2 



state system is of the form: 







where i = 7-i 



and <|> 1= <f*2 are ^bitrary rea l numbers. Unlike the W-H 
transformation, the probability in each state stays the 
same. 

The selective inversion of the phase of the 
amplitude in certain states is a particular case of selec- 
tive rotation that we will need in all of the algorithms in 
this paper. Based on [BBHT], the following is a realiza- 
tion. Assume that there is a binary function/(x) that is 
either or 1 . Given a superposition over states x , it is 
possible to design a quantum circuit that will selectively 
invert the amplitudes in all states where f(x) = 1 . This 
is achieved by appending an ancilla bit, b and consider- 
ing the quantum circuit, as shown in figure below, that 
transforms a state \x, b) into \x, f(x)XOR b) (such a 
circuit exists since, as proved in [Revers], it is possible 
to design a quantum mechanical circuit to evaluate any 
function/(x) that can be evaluated classically.) If the 

bit b is initially placed in a superposition -p(|0) - |1)) , 

this circuit will invert the amplitudes precisely in the 
states for which f(x) = 1 , while leaving amplitudes in 
other states unchanged. 



n bit input: Xj, x 2 . . . x n 



1 



n bit output 
► 



/(*) 



single ancilla bit b 



I 



^ XOR 



single bit output 
► 



4=(io)-ii» 

72 



Figure: The above quantum mechanical circuit 
inverts the amplitudes of precisely those states 
for which the function f(x) is 1. 



3. Amplitude amplification 

Let each point of the domain of f(x) be mapped to a 
state - let t be the target state, i.e. the function f{x) is 
non-zero at the point corresponding to state t. The object 
is to get the system into the f-state. Assume that we have 
at our disposal a unitary transformation U and we start 
with the system in the s-state. If we apply U to s, the 
amplitude of reaching t is U ts , and if we were to 

observe the system at this point, the probability of get- 

2 

ting the right state would be \U t J . It would therefore 



take Q. 



repetitions of this experiment before a 



single success. This section shows how it is possible to 
1 



reach state t in only O 



U. 



steps. This leads to a siz- 



able improvement in the number of steps if \U fs \ « 1 . 

Denote the unitary operation that inverts the ampli- 
tude in a single state x by / . In matrix notation this is 

the diagonal matrix with all diagonal terms equal to 1 , 
except the xx term which is -1 . 

v x denotes the column vector which has all terms 

zero, except for the x' h term which is unity. 

Consider the following unitary operator 1 : 
Q = -I S U l I t U . Note that since U is unitary, U 1 is 
equal to its adjoint, i.e. its conjugate transpose. We first 
show that Q preserves the two dimensional vector space 

spanned by the two vectors: v s and (U l v ( ) (note that 
in the situation of interest, when U ts is small, the two 
vectors are almost orthogonal). 

First consider Qv s . By the definition of Q , this 

-1 T 

is: -I S U I t Uv s . Note that v x v x is an NxN square 
matrix all of whose terms are zero, except the xx term 



1 . According to the notation of this paper (and 
most quantum mechanics texts), the opera- 
tion AB denotes that the sequence of opera- 
tions is B and then A. Q = -I S U l I t U 
implies the following operation sequence: 



First U , then / , then U and then -/ 



T T 

which is 1 . Therefore I f = I - 2v f v f & I s = I - 2v s v s : 
(1) Qv s = -{I-2v/ s )U-\l-2v/ t )Uv s 

= -(I - 2v/ s )U~ l Uv s + 2(1- 2v s v T s )U~\v t vJ)Uv s 



Using the fact that v s v s = 1 , it follows that: 

(2) Qv s = v s + 2(I -2v/ s )U~\v/ t )Uv s 
Simplifying the second term of (2) by the following 

j T — 1 * 

identities: v t Uv t ,= U u and v s U v ( =U ts (the sec- 
ond follows from the fact that U is unitary and so its 
inverse is equal to its adjoint.) 

(3) Qv s = v s (l-4\U ts \ 2 ) + 2UJU-\) 

Next consider the action of the operator Q on the 
vector U l v f . Using the definition of Q (i.e. 

Q = -I U I,U) and carrying out the algebra as in the 
computation of Qv above, this yields: 

(4) Q{ir\)=-i s ir l i t mir\) 

j T — 1 * 

Writing I s as / - 2v v v s and as in (3), v s U v { = U u : 

(5) Q(U~ l v t ) = (U-\ t )-2v/ s (U-\ t ) 

= (U-\)-2U* s v s 
It follows that the operator Q transforms any 
superposition of the vectors v ? & (U l v f ) into another 



superposition of the two vectors, thus preserving the two 
dimensional vector space spanned by the two vectors v s 

& ( U 1 v f ) . As indicated in the figure below, (3) & (5) 
may be written as: 



(6) Q 



U v, 



(l-4|tg )2U tl 
-2U„ 1 



U v, 



It follows as in [BBHT], that if we start with v , then 
after r| repetitions of Q we get the superposition 



a s v s + fl ? t/ v t where a s = cos (2r) I f/ f I ) and 



\a\ = |sin(2Ti|f/J)|. If ti = 



4 U. 



then we get the 



superposition U v t , from this with a single application 

1 



of U we can get v t . Therefore in O 



U 



ts\ 



steps, we 



can start with the s-state and reach the target state t with 
certainty. 

The above derivation easily extends to the case 
when the amplitudes in states, s and t, instead of being 
inverted by I g and I t , are rotated by an arbitrary phase. 

However, the number of operations required to reach / 
will be greater. Given a choice, it would be clearly better 
to use the inversion rather than a different phase rota- 
tion. Also the analysis can be extended to include the 

case where /, is replaced by V l ItV , V is an arbitrary 
unitary matrix. The analysis is the same as before but 
instead of the operation U , we will now have the opera- 
tion VU . 



Q 



(l-4|tg )2U ts 
-2E/,„ 1 



where Q = -I S U I f U 



Figure: The operator Q = -I S U I ( U preserves the 2 dimensional vector space spanned by v s & 

U l v t , it rotates each vector in the space by approximately 2 1 U ts \ radians. Since v s & U l v t 

% -l 
are almost orthogonal, in approximately -j——; applications, it transforms v s into U v f . 



4. Examples of quantum mechanical algorithms 

The interesting feature of the analysis of section 3 is that 
U can be any unitary transformation - furthermore s and 
t can be any two basis states. Clearly, it can be used to 
design algorithms where U is a transformation in a 
quantum computer - this paper gives five such applica- 
tions. The first three are search-type applications, the 
next two are statistical applications. Additional applica- 
tions are presented in [Qntappl]. 

The following general approach is used in each of 
the algorithms. The problem is reduced to one of getting 
the system into a state t. A unitary transform U and the 
initial state s are selected and U fs is calculated. By sec- 
tion 3, it follows that if we start with initial state s and 

carry out o( j^— : J repetitions of the operation 

sequence -I S U l I ( U followed by a single application 

of U , there is an appreciable amplitude of the system 
reaching the f-state. 

The analyses and derivation of the algorithms is 
considerably simpler than comparable analyses and der- 
ivation of existing algorithms. The reason for this is that 
the result of section 3 derives the amplitude amplifica- 
tion based just on U ts and the unitarity of the matrix U, 

i.e. it does not require us to keep track of the amplitudes 
in the remaining states during intermediate steps. 

This power is most evident in the derivation of the 

~J ste P algorithm to estimate the mean to a precision 

£ ((v) in this section). Anyone who has looked at the 
paper [Mean] that derived a, much more complicated, 

step algorithm for the same problem will 
immediately appreciate the usefulness of this approach. 



(i) Exhaustive search starting from the state 

The N = 2 n states to be searched are represented by n 
qubits. In case the starting state s be the state and the 
unitary transformation U is chosen to be W (the W-H 
transformation as discussed in section 2), then t/ for 

any target state t is -J= . This section gives an algorithm 

Jn 

requiring ^(j^ - jj' 0{jN) steps to reach the t- 

state with certainty. Although the derivation is different, 
the end result is the same quantum search algorithm as 
[Search 1] and [Search2]. 

As indicated in section 3, this algorithm starts with 

the with state carries out repeated operations of 
Q = -I s U~ l I t U . With U~ l = U = W and s as the 
state, Q becomes Q = -I-WI.W ; hence the operation 

sequence is: ...(-I-WI f W)(-I-WI t W)(-I-WI t W)... . 

1 ' 

By rearranging parentheses and shifting minus signs, 
this may be seen to be alternating repetitions of 

-WI-W and /,. 
' 

The operation sequence -WI-W is simply the 
inversion about average operation [Search2]. To see 

T 

this, write /- as / - 2v-v- . Therefore for any vector x : 

-WI-Wx = -W(I-2v-v T -)Wx = -x + 2Wv-v T -Wx 
( 0' 

T 

. It is easily seen that Wv-v- Wx is another vector each 
3 o o 

of whose components is the same and equal to A where 

N 

A = — j x- (the average value of all components). 
N 

i = 1 

Therefore the i th component of -WI-Wx is simply: 

(-x ( . + 2 A) . This may be written as A + (A- x { ) , i.e. 

each component is as much above (below) the average 
as it was initially below (above) the average, which is 
precisely the inversion about average [Search2]. 



(ii) Exhaustive search starting from an arbitrary 
basis state 

In case U be chosen as the W-H transform W, the 
matrix element t/ between any pair of states s and t is 

±-|= . Therefore the search algorithm can start with any 

of the N = 2 n basis states as the initial state s and the 
procedure described in section 3, yields an algorithm to 
1 



reach t in O 



U 



ts\ 



, i.e. 0{JN) , iterations. Therefore 



instead of starting with the state, as in (i), we could 
equally well start with any basis state s, and repeatedly 
apply the operation sequence Q = -I WI t W to obtain 

an equally efficient O (JN) algorithm. 

The dynamics is similar to (i); only there is no 
longer the convenient inversion about average interpre- 
tation. The approach of section 3 makes it possible to 
derive this algorithm since we just need to obtain the 
values of t/ for a single application of U , whereas in 

[Searchl][Search2] we had to take the analysis all the 
way to the end based on the inversion about average 
interpretation. 



(iii) Search when an item near the desired state is 
known 

Problem Assume that an n bit word is specified - the 
desired word differs from this in exactly k bits. 
Solution The effect of this constraint is to reduce the 
size of the solution space. One way of making use of 
this constraint, would be to map this to another problem 
which exhaustively searched the reduced space using (i) 
or (ii). However, such a mapping would involve addi- 
tional overhead. This section presents a different 
approach which also carries over to more complicated 
situations as discussed in [Qntappl]. 

Instead of choosing U as the W-H transform, as 
in (i) and (ii), in this algorithm U is tailored to the prob- 
lem under consideration. The starting state s is chosen to 
be the specified word. The operation U consists of the 



transformation 



i - 1 4= 



a 



, applied to each of 



the n qubits (a is a variable parameter yet to be deter- 
mined) - note that if a is 2 , we obtain the W-H trans- 
form of section 2. Calculating t/ it follows that 



n-k 



, this is maximized when a 



is chosen as - ;then log I U ts 



n, n-k k, n-k 
-log n l °g—j— 

2 n 2 k 

The algorithm of section 3 can now be used - as in (i) 
and (ii), this consists of repeating the sequence of opera- 
tions -LU~ l LU, o(v^- 
\\U,, 



times, followed by a single 



application of U . 

The size of the space being searched in this prob- 
n\ 



lem is C, which is equal to 



Using Stirling's 



'* n n-k\k\ 

approximation: logn! = nlogn - n , from this it follows 

that log^C^ = nlog— ^— - fclog— ^— , comparing this to 

Yl — /C Yl — fC 

the number of steps required by the algorithm, we find 
that the number of steps in this algorithm, as in (i) and 
(ii), varies approximately as the square-root of the size 
of the solution space being searched. 



(iv) Estimating the median to a precision e 

Assume that we are given N = 2 n values between 
and 1 denoted by x , x l ...x Nl and a certain thresh- 
old . Let e denote the fractional difference in the 
number of values above and below , i.e. the number of 
N 

values below is — ( 1 - e) . Given the bound |e| < 2e , * ~ 4 Eq 

the task is to find an estimate 8 such that projects the system onto one of its basis states, the prob- 



Also note that U 1 = W l R V 1 = WRW . There- 
fore, by the analysis of section 3, it follows that after 

0\ — | repetitions of Q = -I-U l I 7 U, followed by a 

\e J r o o 

single application of U , the amplitude in the state 
reaches - — . Now if a measurement be made which 



|e||<j. 



If we can solve the above problem in O 



steps, then it is possible to iterate with different thresh 
olds and estimate the median to a precision e in O ( 



steps [Median]. Classically it needs Q.\ -i I steps lo esti- 



ability of getting is 



48, 



It follows from the cen- 



tral limit theorem [Feller] that by repeating this 



experiment O(M) times, it is possible to estimate 



with a precision of -— , and hence |e| with a precision 

Jm 



—— . By appropriately choosing M , it is possible to 

Jm 



mate the median with a precision of e . (All bounds in 
this and the next example are upto polylogarithmic fac- 
tors). 



estimate lei within an error bound of — . 

4 



Solution: This sub-section gives a scheme for the esti- 



mation of e with an error bound of — in O 

4 Ve, 



steps. 



The basic idea is to devise a unitary operation U due to 
which the amplitude in a f-state comes out to be propor- 
tional to the statistic to be estimated. By repeating the 

operation -I S U l I ( U, an appropriate number of times, 

the probability in the f-state is boosted to the point that it 
can be estimated by carrying out the experiment a few 
times. 

n 

Consider an N = 2 state quantum system rep- 
resented by n qubits - associate a value x a with each 

state, S a . Next consider the unitary transform R which 

is a selective inversion operation (as discussed in section 
2): 

R: In case the value x a associated with the state S a is 

smaller than the threshold , invert the amplitude 
in S a . 

Start the system in the state and consider the uni- 
tary operation U = WRW . It is easily seen that after 
U the amplitude of the system being in the state is e . 



(v) Estimating the mean to a precision e 

Assume that we are given N = l" values denoted by 
x , x l ...x N _ l , each x a lies in the range(-0.5, 0.5) - 

the task is to estimate the mean (denoted by |i) to a 
specified precision e . Classically this would take 

^(""2) steps. The best known quantum mechanical 
algorithm takes (9^— j-^j steps [Mean]. This section 

presents a simple ^f~J step quantum mechanical algo- 
rithm (upto polylogarithmic factors). 



states. We need the following 4 unitary operations in 
this algorithm (these are denoted by My ,M 2 , Ry and 

Wy .) It may be verified that these are indeed valid uni- 
tary operations and, with the above encoding, they 
require operations on only a few qubits at a time. 

Mj If in state S : go to Q with an amplitude of 
1 



; , stay in S with an amplitude 
73 . 



of 



if in state Q: go to Sq with an amplitude of 
- , stay in Q with an amplitude 



Start with a relatively large e which is chosen so 
that < e . Carry out the following loop (i)...(iii): 

(i) Estimate (0. with a precision of e () , i.e. find an 

e 

estimate \i e such that - |i| < — . This step is 

described in the "Main algorithm" below. 

(ii) Shift each of the numbers by the newly estimated 
mean, i.e. x a = x a -\l g . 

e 

(iii) If £ () > e , replace e () by — and go to (i). 

The mean may be estimated as the sum of the estimated 
\l in each iteration of step (i) of the loop. 



of-^; 
2 

if in any 

other state: stay in the same state. 
M 2 If in state S : go to Q with an amplitude of 

, stay in S with an ampli- 

tude of 4= ; 

if in state Q: go to S with an amplitude of 
1 



72 



, stay in Q with an ampli- 



Main algorithm The heart of the above algorithm is 
step (i) of the above loop, i.e. estimating \i with a preci- 
sion of e when < e () . Consider a (2 n+ 1 + 1) state 

quantum mechanical system represented by (n + 2) 
qubits with the following encoding for the states: 

So Sj S NA 

(first two bits are 00, next n bits indicate the S a state) 



R Rj 



R 



N-l 



(first two bits are 01, next n bits indicate the R a state) 

Q 

(first bit is a 1, the next (n + 1) bits are zero). 
The S a states are the computational states, for each S a 
state there is an R a state, also there is a single Q state. 
Each state is encoded by (n+2) qubits as shown above. 



Associate the value x a with each of the R a and S a 



tude of — ; 

72 

if in any 

other state: stay in the same state. 
R x If in state S a : go to R a with an amplitude of 



2 4 2 

- - -x a , stay in S a with an 



( 1 2ix a 
amplitude of — + — — , , 

V73 73 J 

if in state R a : go to S a with an amplitude 



ofJ^-^ a 2 , stay in R a with 

( 1 2ix a 
an amplitude + 

V 73 73 

if in any 

other state: stay in the same state. 



W x If in state S a : carry out the W-H transform on 
the S a states; 

if in any 

other state: stay in the same state. 

It is easily shown that if we start with the Sq state and 
carry out the operation U , where U is defined as: 
U = M l W l R l W l M 2 , then the amplitude in the S () 

state will be ^ . As in the median estimation, (iv), U 1 

is the product of the inverses of the same matrices as U , 
but in the opposite order. By the analysis of section 3, it 

follows that: By o( — ] repetitions of Q where 



Q = -I S JJ I s U followed by a single application of 



U , the amplitude in the S state becomes . Now if 

V2eo 

a measurement is made which projects the system onto 
one of its basis states, the probability of getting the Sq 

state is \( — 1 . It can be shown by the central limit the- 



2le 



orem that by repeating this entire experiment O(M) 
times, it is possible to estimate — with a precision of 



4m 



By an appropriate choice of M , it is possible to 



estimate (J. within an error bound of — . 

Note that in the above analysis we only needed to 
design the unitary transform U to appropriately adjust 
the one matrix element that we cared about - just the 
condition that it was unitary took care of everything 
else. 



5. General quantum mechanical algorithms 

The framework described in this paper can be used to 
enhance the results of any quantum mechanical algo- 
rithm. Assume there is a quantum mechanical algorithm 
Q due to which there is a finite amplitude Q u for transi- 
tions from the starting state s to the target state t. The 

2 

probability of being in the state t is hence \Q - it will 



therefore take Q. 



\Q ts \ 



repetitions of Q to get a single 



observation of state t. Since the quantum mechanical 
algorithm Q is a sequence of r| elementary unitary 
operations : Q Y Q 2 . ■ . , it is itself a unitary transforma- 
tion. Also, Q 1 = Q^Q^ _j ---Qi ■ The inverse of any 

elementary unitary operation on a small number of 
qubits is another elementary unitary operation on the 
same qubits and can hence be synthesized. Applying the 
framework of section 3, it follows that by starting with 
the system in the s state and repeating the sequence of 



operations: 



times followed by a 



single application of Q , it is possible to reach the f-state 
with certainty. 

6. Conclusion 

Designing a useful quantum computer has been a daunt- 
ing task for two reasons. First, because the physics to 
implement this is different from what most known 
devices use and so it is not clear what its structure 
should be like, The second reason is that once such a 
computer is built, few applications for this are known 
where it will have a clear advantage over existing com- 
puters. This paper has given a general framework for the 
synthesis of a category of algorithms where the quantum 
computer would have an advantage. The main concept is 
that if there is a finite amplitude for a certain transition, 
this can amplified by a factor of M in M iterations of 
the algorithm. Several applications were presented. It is 
expected that this formalism will also be useful in the 
physical design of quantum computers, since it demon- 
strates that quantum mechanical algorithms can be 
implemented through general properties of unitary 
transformations and can thus adapt to available technol- 
ogy. 
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